Multiple regimes of diffusion 
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We consider the diffusion of independent particles experiencing random accelerations by a space- 
and time-dependent force as well as viscous damping. This model can exhibit several asymptotic 
behaviours, depending upon the limiting cases which are considered, some of which have been 
discussed in earlier work. Here we explore the full space of dimensionless parameters, and introduce 
an 'asymptotic phase diagram' which delineates the limiting regimes. 
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I. INTRODUCTION 

The position x of a particle subjected to a force which 
fluctuates randomly in time t might be expected to un- 
dergo diffusion, in the sense that 
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(1) 



for some diffusion coefficient V^, provided that there is 
some damping mechanism preventing the particle from 
being accelerated to arbitrarily high velocities. In this 
paper we determine the diffusion constant for the sim- 
plest model for this process in one spatial dimension, in 
which the equation of motion of the particle is 
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Here x and p are particle position and momentum, re- 
spectively, m is the mass, and 7 is the rate at which 
the particle momentum is damped due to viscous drag. 
Time derivatives are denoted by dots. Further, the ran- 
dom forcing f{x,t) is modeled by a Gaussian random 
function with zero mean and with correlation function 
C{x,t), characterised by a correlation length ^, correla- 
tion time T, and of typical size a: 



{f{x,t))^0, {f{x,t)f{x',t'))^C{a 



-t') (3) 



(angular brackets denote averages). The dynamics of the 
model ([2]), ([3]) is determined by five dimensional param- 
eters: cr, r, ^, the mass to, and 7. Out of these one 
can form two independent dimensionless parameters: a 
dimensionless force x = crT'^/{m^) and a dimensionless 
damping u = 7T. In the following we explore the full 
space of dimensionless parameters x ^-nd lo. 

There is no exact expression for the diffusion constant 
for this simple model. However, asymptotic expressions 
with different regions of validity are known, depending on 
which ratios of dimensionless parameters, uj and Xi ^P" 
proach zero. We show that there are surprisingly many 
different asymptotic regimes, which are summarised in 
an 'asymptotic phase diagram'. Fig. 1. The axes in this 
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FIG. 1: Phase diagram (schematic) for the model ([2]|4|) sum- 
marising the different dynamical behaviours of Eq. ((2} de- 
scribed in sections IIIIllVl I Ornstein-Uhlenbeck, II gener- 
alised Ornstein-Uhlenbeck, Illb overdamped minimum track- 
ing. Ilia underdamped minimum tracking. The Ornstein- 
Uhlenbeck regime is divided into three regions: overdamped 
advection (la), underdamped advection (lb), and under- 
damped inertial dynamics (Ic). 



diagram are logarithms of two independent dimensionless 
parameters of the model, and a ray from the origin with 
slope v represents a limiting process where x approaches 
or 00 with Lu ^ x'^ ■ The phase lines do not indicate 
sharp transitions, but rather the boundaries between the 
domains of validity of six different asymptotic regimes as 
the limit is taken. Some of the regimes are well under- 
stood, but others are either new or have only been studied 
recently by the authors of this paper. It is remarkable 
that the phase diagram of such a fundamental model for 
diffusion processes has not been completely characterised 
before now. In all cases the long-time dynamics is dif- 
fusive, but in some of the regimes the stationary dis- 
tribution of momentum may be strongly non-Gaussian, 
and the short-time behaviour may exhibit anomalous dif- 
fusion. The diffusion constant depends in different 
ways on the microscopic parameters in different regimes, 
as illustrated in Fig. O 
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FIG. 2: Illustration of multiple regimes of diffusion. Evalua- 
tion of the diffusion constant for constant x as a function of 
uj: for X = 50 (circles) and x = 1250 (squares). Also shown is 
the expected behaviour in the advective limit, la, (blue) and 
the expected behaviour in regime II, T>x ~ ^^7''" X^^^'-^~^^^ 
(red). Finally, the estimate (|42p for "D^ in the minimum- 
tracking regime is shown (green line). 



The three-dimensional version of the model defined by 
Eqs. (HI) and ([3]) arises naturally in the study of small par- 
ticles suspended in a randomly moving fluid, for which 
motion relative to the fluid is determined by viscous drag. 
In that context the random force is replaced by a ran- 
dom vector field, which would usually be chosen to be 
solenoidal, to represent an incompressible flow. This 
three-dimensional system has been extensively studied 
in certain limits. A significant early contribution is due 
to Maxey [l|, who analysed the clustering of particles 
suspended in a turbulent fluid (referred to as 'preferen- 
tial concentration'). Ref. provides an overview of the 
literature on this problem, and describe recent progress. 

This present paper is the first to explore the full range 
of regimes which are possible in limiting cases of the 
model, some of which are not realised in fluid-dynamical 
applications. We remark that different choices of dimen- 
sionless parameters are used in some other papers: much 
of the fluid dynamics literature uses the Stokes number 
St = l/cj as a measure of the damping, and the Kubo 
number Ku = x/w as a measure of the time scale of 
fluctuations of the velocity field. 

The model also exhibits an interesting effect which in- 
volves a phase transition in the conventional sense. De- 
pending upon the dimensionless parameters of the model, 
particles with different initial conditions experiencing the 
same realisation of the random force approach the same 
trajectory with probability unity. This 'path coalescence' 
effect and the 'path-coalescence transition' where it dis- 
appears, were noted by Deutsch who appears to have 
been the first to consider this model systematically. In 
this paper we also describe the full phase line for the 
path-coalescence transition, extending results of Q . The 
critical line for the path-coalescence transition in the 
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FIG. 3: Phase line of the path-coalescence transition. Results 
of computer simulations of (OH for the phase line of the path- 
coalescence transition. Also shown is the theoretical result 
valid for small values of x, red line. 



plane is shown in Fig. [3) 

The numerical simulations of Eqs. ([2]), Q described in 
this paper were performed with the following choice of 
correlation function: 



C(x,t) 



tj^ exp [ 



(4) 



The detailed choice of the correlation function of the force 
is not significant, but in regime II, one aspect of the ran- 
dom forcing can make a qualitative difference. For any 
choice of the random force, there is a corresponding po- 
tential, satisfying —dV{x^t)/dx = f{x,t). For a generic 
choice of correlation function, the one-dimensional poten- 
tial V{x,t) corresponding to the force f{x,t) performs 
a random walk exhibiting increasing fluctuations as 
increases. We also consider cases where the particle dy- 
namics is different if the potential V{x, t) is a stationary 
random process. 

The different regimes are illustrated in Fig. 3] by nu- 
merical simulations of the Eqs. ([2J3]). Shown are the tra- 
jectories x{t) of several particles for a given realisation 
f{x,t) of the forcing. 



II. SUMMARY AND PHYSICAL DESCRIPTION 
OF THE REGIMES 

Before describing our results in detail we discuss the 
physics of the parameter regimes of our model. We also 
mention connections with other work on the dynamics of 
randomly forced particles, where some of the regimes of 
our model have been studied. 

One limiting case is where the particles are advected 
by a random velocity field: there is extensive literature 
on this problem and the closely related model of passive 
scalars [5|. The advective case corresponds to regime la 
in our model. 
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FIG. 4: Trajectories x{t) as a function of t/r for 20 particles 
with initial condition p(0) = and a:(0) randomly chosen 
in [0,1]. The first two pictures (in regimes la and Ic) are 
similar to those obtained by Deutsch jsj, and to Fig. 1 in [j. 
The trajectories in region lb are similar to those in la (not 
shown). The general difference in the dynamics of regimes 
Ilia and Illb is that there are oscillations in regime Ilia, due 
to smaller damping. 



Our model ^ , ^ reduces to the well-known Ornstein- 
Uhlenbeck process [gj] when the position dependence of 
the force can be neglected. This is the case when the 
forcing is sufficiently weak so that the particle position 
changes negligibly within correlation time r: this condi- 
tion will be made more precise below. The Ornstein- 
Uhlenbeck process is discussed in standard textbooks 
(see for example Q)- The regimes Ib,c can both be 
analysed by treating the evolution of the momentum as 
an Ornstein-Uhlenbeck process. The difference between 
the two regimes is that regime lb exhibits path coales- 
cence whereas regime Ic does not. Despite the fact that 
regime la describes an overdamped process, the formula 
for the diffusion constant is the same as for the Ornstein- 
Uhlenbeck process. For this reason regimes la and Ib,c 
are treated together in Sec. lIIIi and are referred to as the 
Ornstein-Uhlenbeck regimes. 

Our model is also related to stochastic or 'Fermi' ac- 
celeration of classical particles by random forces, which 



is used to model the production of cosmic rays [8|. In 
these studies the damping term (proportional to 7) is 
not included in the equation of motion, and the particle 
is accelerated to arbitrarily high energies. The treatment 
of the random forcing term in the case where the particle 
is rapidly moving, first considered in Q, is used in our 
consideration of regime II. Without damping the model 
exhibits anomalous diffusion 10, JTj. For the case where 
the damping term is included, a new dynamical regime 
was identified in [l^, [l^ with a non-Maxwellian veloc- 
ity distributions (as well as anomalous diffusion at short 
times, before the damping term starts to limit the accel- 
eration). This is regime II in the phase diagram Fig. [TJ 
we call it the 'generalised Ornstein-Uhlenbeck regime'. It 
is discussed in Sec. ITVl 

In the case where both the damping and the force are 
strong, the particle follows a local minimum of the po- 
tential V{x^t). This is regime III of the phase diagram. 
This 'minimum-tracking' regime has not been considered 
in detail in earlier work. It is discussed in Sec. IVl below. 

In Sec. lVII we briefly describe how results for regime II 
differ for more general types of forcing, such as the case 
where V{x,t) has stationary statistics. 

Fig. [3] shows numerical results on the path-coalescence 
transition, for the choice of correlation function For 
small values of x, the phase boundary is in precise agree- 
ment with an asymptotic theory discussed in Q, which 
shows that the transition line is determined by the con- 
dition wx"^/^ const, in the limit as x ^ 0. The data 
for large x are consistent with the hypothesis that the 
phase line approaches oj — const, as x ~> 00, but we have 
no compelling argument to support this. 

Finally, we comment on the physically accessible range 
of dimensionless parameters. This depends on the na- 
ture of the forcing. In the case of a particle suspended 
in a turbulent fluid flow with velocity field u{x,t), the 
random forcing is due to viscous drag, and we write 
f(x,t) = m'yu{x,t). In this case disturbances in the 
fluid velocity field m(x, t) are transported by u{x, t) itself. 
This implies that the Kubo number, Ku = ur/^ cannot 
be large (equivalently, x cannot be large compared to uj) 
if the random forcing is due to a fluid flow. In other 
cases, such as forcing by random electromagnetic fields, 
the entire phase diagram may be accessible. 



III. ORNSTEIN-UHLENBECK REGIME 

We term the regimes Ia,b,c in the phase-diagram Fig.[T] 
the 'Ornstein-Uhlenbeck regimes'. In these regimes, the 
particles move so slowly that the distance travelled dur- 
ing one correlation time r of the random force f{x,t) is 
much smaller than its correlation length ^. Thus changes 
in the spatial argument do not contribute significantly 
to the decorrelation of f{x,t) and one may approximate 
f{x{t),t) w f{x{0),t) for times t of the order of or less 
than the correlation time r. 

Regime I is divided into one overdamped regime, 
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regime la {lo ^ 1) and two underdamped regimes lb, 
c (w < 1). 

In the overdamped case, the acceleration term in eq. 
(m is negligible, and consequently the particles are ad- 
vected by the random force, that is, x ~ f (x , t) / (rwy) . 
The two regimes lb and Ic are distinguished by different 
behaviours of nearby particles, as can be seen in Fig. 31 
In regime lb, initially separate but nearby particle trajec- 
tories approach each other (path-coalescence regime). In 
regime Ic, by contrast, initially close particle trajectories 
do not coalesce. 

In the remainder of this section we first briefly describe 
the single-particle dynamics (diffusion), and then sum- 
marise what is known about the path-coalescence transi- 
tion. 

To describe spatial diffusion in the overdamped regime, 
one integrates the advective equation of motion x ~ 
f {x , t) / (rwy) . The change in position 6x during a short 
time interval St ^ t is 



Sx 



2 rt+St 
Jt 



(5) 



In regime I, one may approximate f{x{ti),ti) w 
f{x{t), ti), as pointed out above. In this regime, the fluc- 
tuations of / at a given point x is indistinguishable from 
the fluctuations of / along a particle trajectory. In this 
case it is straightforward to determine the fluctuations 
of Sx: since the force is assumed to have vanishing mean 
([3|), one has (Sx) = 0. The variance of Sx is determined 
by making use of the fact that (/(t)/(0)) is small unless 
\t\ < T. Evaluating (Sx^) one finds the standard result 



{Sx^ 



where Dq is given by 



2DoSt 
(7717)^ 



dt{f{t)fm . 



(6) 



(7) 



The position x{t) at time t — NSt of a particle after 
N microscopic steps is x{t) — x{0) = '^i^iSx'-'^\ where 
fe^^^ is the increment at the time step number i. For 
the spatial diffusion constant one obtains in the usual 
fashion: 



— lim 

t — >oc 



{{x{t)-x{o)r) 



2t 



(8) 
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2StN 



where the increments fe'*-* and Sx^^'> are statistically in- 
dependent when i ^ j and St ^ r. 

Consider now the underdamped regimes lb, c. The 
displacements Sp of momentum (for a short time interval 
(5t) obey 



Sp = —^pSt + Sw , 



(9) 



with 



Sw 



t+st 



dti f{x{ti),ti 



(10) 



In regime I, the force fluctuates sufficiently rapidly com- 
pared to the time scale on which the momentum relaxes 
(lu ^ 1) and the change Sp of momentum during one cor- 
relation time of the force is small compared to the typi- 
cal value of p. Consequently ([5]) is a Langevin equation 
01, describing the standard Ornstein-Uhlenbeck process 
where Sw is Gaussian distributed with 



(Sw) = , (Sw^) = 2DoSt . 



(11) 



From the corresponding Fokker-Planck equation for the 
distribution P{p, t) of momentum p at time t [7| 



dP d , d . 



(12) 



one deduces that the steady-state distribution of momen- 
tum is Gaussian, P{p) oc exp(— 7p^/2£)o)- 

Also in this underdamped regime the particles diffuse. 
The Fokker-Planck equation (jl2p allows to determine the 
correlation function of momentum in the steady state. 
The result is [7]: 



steady 



^0 



state - ^ exp(-7|t2-tl|). (13) 



This result in turn allows to calculate the spatial diffusion 
constant: 



i.„. 'fm 

t— ►oo 2t 

1 



(14) 



lim 
t^oo 2im^ 



t rt 
dt 



dt2 {p(tl)p{t2)) 

steady f 



(771,7)^ 



By comparing the results (|8]) and p4)) one finds that the 
spatial diffusion constant is the same in the over- and the 
underdamped limits. 



Do 
(777.7)^ 



(15) 



Fig. [2] shows results of numerical simulations for the dif- 
fusion constant of the model ([2111]). In the Ornstein- 
Uhlenbeck regime (I) the simulations agree well with 
Eq. (Uni). 

We now briefly summarise what is known about the 
path-coalescence transition which distinguishes regime lb 
from Ic. In the path-coalescing phase (regime lb), par- 
ticle trajectories governed by the equation of motion ([2]) 
coalesce, whereas in regime Ic initially close particle tra- 
jectories separate almost surely (see Fig. S]). As was ar- 
gued in 4], the maximal Lyapunov exponent A serves as 
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an 'order parameter' for the phase transition. The ex- 
ponent describes the rate of change of an infinitesimal 
separation between two trajectories 



A = Um t ^ logp 



Sxf 
6x0 



(16) 



Here Sxq is the initial separation of two infinitesimally 
close trajectories, and 6xt is their separation at time t. 

In regime lb, the Lyapunov exponent is negative while 
it is positive in regime Ic. The condition for the phase 
transition is thus A = 0. In regime I, the Lyapunov expo- 
nent can be calculated exactly [1| . Expressed in terms of 
the dimensionless parameters x and ui, the phase transi- 
tion is found to occur at 



u;X 



-2/3 _ 



const. 



(17) 



Fig.[3]shows results of numerical simulations for the locus 
of the path-coalescence transition in the x-w-plane for the 
model (I2E1). In regime I, the transition line is given by a 
line of slope 2/3, as expected from pT|) . 

To conclude this section, we briefly discuss the condi- 
tions delineating regime I in Fig. [T] First, in the over- 
damped limit, it was assumed that the distance covered 
in time r is smaller than ^. Estimating the advective 
velocity by (7/(7717) have the condition ar / {m'y) <C 
Consequently the line distinguishing regimes la and Illb 
is given by the condition that w/x is of order unity. Sec- 
ond, in the underdamped limit, w <^\, the boundary be- 
tween regimes Ic and II is parameterised by the condition 
that oj/x^ is of order unity. This condition is derived in 
the next section (Eq. where regime II is discussed. 



IV. 



GENERALISED ORNSTEIN-UHLENBECK 
REGIME 



We term the regime II in the phase diagram Fig. [T] the 
'generalised Ornstein-Uhlenbeck regime'. This regime is 
defined by strong stochastic forcing: the particles move 
fast and their momenta can be much larger than po = 
777^/r. This means that the particle may travel many 
correlation lengths ^ during one correlation time t and 
the stochastic force acting on the particle may decorrelate 
quicker than r. Thus the 'effective correlation time' of 
the force ^m/ \J (p^) can be much smaller than r. It is also 
assumed that the dynamics of regime II is underdamped, 
i.e. w <C 1. 

A Fokker-Planck description is adequate in this regime. 
Define the increment 5w of the force for a small time 
interval 5t 



Sw — 



t+st 



dh f{x{ti),h) 



(18) 



In regime I, the time dependence of x{ti) could be ne- 
glected, but this approximation is no longer valid when 



the forcing is strong. Instead, one integrates the equation 
of motion ^ to obtain x{t) — x{0) + 5x, where 

Sx=- /" Ve-T*' [p(0) + f'dt"e-''"f{x{t"),t")\. (19) 
mJo \- Jo -I 

For times smaller than the small time interval 6t, Sx is 
small and one may expand the force around 6x — 0. To 
lowest order in St one obtains 



(Sw) 



dh{j-{x{{))M)5x 



x(^(x(0),ii)/(x(t"),t") 



**dti f'dt' f dt"e-'(*'-*") (20) 
Jo Jo 

'A 

dx 



where we have used that {f{x{0),t)) = 0. Since it is 
assumed that uj = jt 1, the exponent j{t' — t") k, 
and the major contribution to the integrals from the force 
correlation is for \ti ~ t"\ < t. We get 



^ ' 2m \dx 



{0,0) fipt/m,t) 



(21) 



The variance of the displacements becomes lowest 
order in 5t 



{6w^) = 



St pSt 

dh dt2{f{pti/m,h)f{pt2/m,t2)) 



Jo 

00 



St / dtC{pt/m,t) 



(22) 



where C is the correlation function of the force ^ . 

Using that the change Sp of momentum during a short 
time period is 5p = —^pSt + Sw together with the first 
two moments of Sw (Eqs. (|2ip and ([H])), a Fokker-Planck 
equation is obtained using the standard procedure Q: 



Here the drift- and diffusion-coefficients are 

d 



(23) 



v{p) 
Dip) 



lim M . 

st^o St 

.5*^0 2St 



'IP + 'q^D{p) 



dtC{pt/m,t), (24) 



where it was used that {Sw) = StdpD{p). The above ex- 
pression for the diffusion constant was earlier obtained 
in [§], and has also been used in [l^ and [ill- Note 
that when \p\ <C po, D{p) « Dq, which corresponds 
to the Fokker-Planck equation of the standard Ornstein- 
Uhlenbeck process discussed in Sec. IIIIl 

On the other hand, when \p\ 3> po we approximate 



D{p). 



DiPo 
\P\ 



/OO 
dXC{X,0). (25) 
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FIG. 5: Non- Gaussian distribution of momentum in phase II 
(c.f. Fig. [1]). Shown are results of numerical simulations of 
dH) for X = 50 and u = 0.01 compared with eq. (PS)). 



For the correlation function ^ we obtain D{p) = 
Do/{l+p^/piy^^, that is Di = Do. 

If the force is the gradient of a potential V{x, t) with 
continuous derivatives, Di vanishes and D{p) oc \p\^'^ 
provided V{x,t) is sufficiently diffcrentiable. Further, if 
the correlation function exhibits a cusp at t = (an 
example is discussed by [13]), we find D{p) oc \p\^^. In 
general we write 



D{p) = D^ 



(26) 



with C ^ 0. In this paper we mainly discuss the case 
= 1 which is a generic case for a random force [realised 
by the correlation function (|3|)]. But in Sec. lVII we briefly 
mention what is known for other force models giving rise 

to C 7^ 1. 

The steady-state solution P{p) of is given by 



Pip) = C exp ( - 7 



dp' 



P' 



Dip'] 



(27) 



where C is chosen to normalise the distribution. When 
Dip) « Do, the function Pip) is approximately Gaussian, 
which corresponds to the standard Ornstein-Uhlenbeck 
regime (regimes Ib,c). When Dip) = Dipo/\p\, however, 
the distribution Pip) is non-Gaussian (regime II): [12] 



Pip) - exp - 



j\pI_ 

3poDi 



(28) 



This result is compared to results of numerical simula- 
tions of Eqs. (04]) in Fig. El 

The generalised Ornstein-Uhlenbeck regime was stud- 
ied in Ref. [13, see also Ref. [11]. The Fokker-Planck 
equation l[23)) was solved as an eigenvalue problem and 
the propagator to reach momentum p at time t given 
the initial momentum was found. From this propagator, 
the momentum correlation function at equilibrium was 





I 10 1001000 

t t 

FIG. 6: Shows {x^)t as a function oi t; si X = 
50 and Lo = 0.025,0.075,0.25,0.5,1,2.5,50,200,500 and 
800, from top to bottom; b x = 1250 and uj — 
0.125, 1, 0.5, 2.5, 5, 50, 100, 250, 1250, 2500 and 5000, from top 
to bottom. Solid lines are fits (as judged by the eye) to the 
diffusion law ([l}. The corresponding values of Dx are shown 
in Fig. H 



calculated 



{pit')pin)e,. 



r(4/3) fpoDi 



3i/3r(5/3) V 7 
V3' 3' 3' 



2/3 



'2j\t'~t"\ 



37l*'-t"l 



(29) 



where i<2i is a hypergeometric function 114). 

From the momentum correlation function (|29p it is pos- 
sible to calculate the spatial diffusion constant in regime 
11. The result is 421: 



iPoD^ni-^' 

771,2^5/3 2 



F /i i 2 5 5 ^ 
r(2/3)2 'H3'3'3'3'3' /' ^ ^ 



where F32 is a hypergeometric function. Since Di 
the spatial diffusion constant "Dx scales as 



(31) 



as opposed to (llSp . Fig.[2]shows results of numerical sim- 
ulations of the diffusion constant for the model in 
this regime, in good agreement with eq. (|30p . The diffu- 
sion constant 'Dx was numerically determined by estimat- 
ing (x2)t/(2t), according to Eq. ([1]). The corresponding 
data are shown in Fig.[6l The dynamics in the generalised 
Ornstein-Uhlenbeck regime exhibits anomalous diffusion 
at short times , Ts^ . Results of numerical experiments 
exhibiting anomalous diffusion are given in [l2l ]. 

To conclude this section we discuss the conditions un- 
der which the results described above are applicable. 
First we discuss the lines in Fig. [1] defining the limits 
of regime II. The discussion above assumes that the dy- 
namics in regime II are described by a Fokker-Planck 
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equation, (|23|) . and we must also consider the conditions 
under which this equation is apphcable. 

Regime II is defined by the condition that the motion 
is underdamped, w 1, and by the condition that the 
correlation time along the particle trajectory is smaller 
than the correlation time for a static particle. The latter 
condition defines the transition between regimes I and 
II in the phase diagram Fig. [TJ This transition occurs 
when ^Pq/Dq is of order unity. This fact is most easily 
seen by determining the steady-state distribution P{p) of 
momentum for the particular form Q of the correlation 

function: P{p) = Ccxp ( - ^ [(1 + p^/plf/^ - l]). In 
terms of the dimensionless parameters lu and x the above 
condition becomes: 

IPq im'^S.'^ , .„„^ 

— ^ = — = const . (32) 

Thus the lines distinguishing the boundaries of regime II 
in the phase diagram in Fig. [1] are 

uj^Om and ^=0(1). (33) 

X 

Consider finally the conditions of validity of ((23)) in 
regime II. For a Fokker-Planck description of a stochastic 
process to be valid, two necessary conditions must hold 

• Amplitude condition. The random jumps of the 
stochastic variable must be much smaller than its 
typical size. Therefore we must require that the 
change of momentum Ap during a correlation time 
of the forcing is much smaller than (p^). This 
condition can be written in the form 

- Xi—j < 1- (34) 

p p p {p'^) ^X 

Here we used the fact that the correl ation time 
of the force is of the order of m^/-^ (p^), and 
that in the steady state ^ {PoDihY^^ ^ 

Pq{x^ /ujY/'^ . The amplitude condition is thus ful- 
filled provided uP' jx^ 1- 

• Frequency condition. The stochastic forcing must 
fluctuate more rapidly than the stochastic variable. 
The c orrela tion time of the force is of the order of 

/ \/ kp^) J and the relaxation time of p is of the 
order of 7""'^. We must therefore require that 

^^^^ — «!■ 35 

The frequency condition thus amounts to the same 
condition as above. 

In our model, the amplitude and frequency conditions 
may not be sufficient to ensure that ([23]) is valid. The 
reason is that the fluctuations experienced by the parti- 
cle may be influenced by the random force altering the 
trajectory of the particle, so that the trajectory does not 




X 



FIG. 7: Disappearance of a minimum in V{x,t) (black lines) 
and corresponding coalescence of particle trajectories (red 
lines, schematic). 

explore the random force field ergodically. In addition to 
the two conditions above the following self-consistency 
condition must also be satisfied: 

• Self- consistency condition. It is required that in the 
steady state the particle moves sufficiently rapidly 
so that it is not captured by 'valleys' in the poten- 
tial corresponding to f{x,t): 

(p2)/(2TO) >^(j. (36) 

This condition too corresponds to uj'^ /x ^ 1- 

Within the boundaries of region II we have w <C 1 
and uj/x^ 1- In this limit the condition ui'^/x <C 1 is 
always satisfied, so that the Fokker-Planck equation is in 
fact always applicable in regime II for the type of random 
force model with C = 1- In Sec. IVII we shall see that this 
may not be true for other values of C- 

V. MINIMUM-TRACKING REGIME 

We term the regimes III in the phase diagram Fig. [1] 
as the 'minimum-tracking regimes'. These regimes are 
defined by large damping a; 3> 1 and strong stochastic 
forcing x/i^ 3> 1. If the force had been weak, we would 
be in the advective regime la. When the force becomes 
large enough, the particle may get stuck in minima of the 
potential V{x,t) = — J^dx'f{x',t) of the force. 

The minimum-tracking regimes are divided into two 
distinct dynamical regimes: under- and overdamped 
minimum-tracking. In the underdamped regime Ilia, 
the particle can oscillate around the potential minimum, 
whereas in the overdamped regime Illb such oscillations 
are quickly damped out (Fig. [5]) . 

In regime Illb the diffusion constant is estimated as 
follows. Typically the particles are stuck in minima. As 
these minima randomly disappear and appear, the par- 
ticle trajectories jump and may coalesce (Fig. [7]). We 
argue that the diffusion constant is given in terms of the 
typical size of the valleys (^) and the rate at which they 
disappear. This rate can be estimated using a general- 
isation of the method discussed in [ll, [l^ for counting 
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the zeroes of a random function F{y). According to p^ . 
the density g of zeroes of the random function F[y) is 
given by 



Q= (<5(F)|detG| 



(37) 



with Gai3 = dFa/dyp. 

Valleys disappear at inflection points of V{x,t), and 
we thus need to count the joint zeroes of / and dxf = /': 



W)'5(/')|dctG|) 



with 



G 



■I dt 

f" §IL 

J dt 



(38) 



(39) 



We need to calculate the average with respect to the set 
of Gaussian random variables = {/,/',/", 9*/, i9t/'}- 
The expectations of all five random variables vanish, and 
their covariances can be expressed in terms of deriva- 
tives of the correlation function C{x,t). For the special 
case (HI we obtain the covariance matrix 



/ a2 -aye 

' aye 

-aye 3aye 



\ 



cryier') / 



Performing the Gaussian average in p8|) . we obtain 

V2 1 



(40) 



(41) 



Note that this result is independent of a. For the choice 
of correlation function ([4]), the density of minima of 
V{x, t) for any given time is {2tt£,)~^. Their typical sepa- 
ration is thus 2tt£^. Assuming that particles jump by 27r^ 
when their minima disappear, we arrive at the following 
estimate for the diffusion constant 



2 TT^ r T 



(42) 



This expression is in reasonable agreement with the re- 
sults of numerical simulations shown in Fig. O 

Now consider the conditions delineating regime Illb. 
As u is increased, regime la is entered from Illb. The 
corresponding crossover line was determined in Sec. IIIIl 
it is given by the condition that w/x is of order unity. On 
the other hand, when u is decreased, regime Ilia is en- 
tered, the underdamped minimum-tracking regime. The 
corresponding cross-over line is given by the condition 



7^valloy 



1 



(43) 



where Wvaiicy = \/ c/(^?7i) is the typical frequency of os- 
cillation in the minimum. The above condition corre- 
sponds to 



Lo'^/x = const. . 



(44) 



Both minimum-tracking regimes are on the path- 
coalescence side of the phase diagram Fig. |3l correspond- 
ing to negative Lyapunov exponent. Calculating this ex- 
ponent is complicated by the fact that the minimum- 
tracking regimes are non-ergodic. However, in regime 
Ilia the Lyapunov exponent can be estimated as follows. 

Making use of the fact that the particles are con- 
strained to follow potential minima of the random force 
f{x,t), expand the potential around a specific mini- 
mum at xq at a given time Iq, V{x,to) ~ V{xo,to) + 
lV"{xo,to)(x - xo)2. This gives / = -d^V « -C{x - 
Xq), where C = V"{xo, to) > 0. In the vicinity of xq, the 
equation of motion ^ becomes 



mx — p, 



P = -1P-C{X - Xq) 



Linearising for small separations 5x and 5p gives 



5x — Sp/m, 



Sp = —jSp — C5x . 



(45) 



(46) 



The dynamics of the separations is determined by the 
eigenvalues of the matrix 



1/m 

-C -7 



(47) 



The eigenvalues of A are 




(48) 



In regime Ilia, the discriminant in (|48p is always nega- 
tive (this follows from (|^^ and the fact that C ^ aj^ 
which gives C/m'^'^ ~ 3> 1). In this regime, the 

real part of the eigenvalue, which describes the stability 
of the minimum, therefore always takes the value —7/2. 
Because the real part of the eigenvalue is constant, we 
can conclude that the Lyapunov exponent in regime Illa 
is A = -7/2. 

Finally, as uj is reduced further, eventually region II 
is entered when the Lyapunov exponent turns positive. 
The phase line in Fig. [3] appears to be consistent with 
the path-coalescence transition occurring at u; = const.. 
However, we have not been able to find an argument sup- 
porting this observation. Moreover, we have performed 
numerical simulations of the path-coalescence transition 
at large values of x for random forces with correlations 
different from ([4]) and have found that the locus of the 
path-coalescence transition sensitively depends on the 
nature of the fiuctuations of the force. This observation 
is in stark contrast to the boundary between regimes Ib,c 
where the scaling of the zero of the Lyapunov exponent is 
universal. This latter universality is illustrated by recent 
analytical results of Falkovich et al. (corresponding 
to regimes Ib,c) for a model where the force gradient is 
taken to be a piece-wise constant function (a so-called 
'telegraph process'). 
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VI. A SPECIAL CASE: GENERALISED 
ORNSTEIN-UHLENBECK REGIME 



Self- consistency condition. For general values of C, 
the self-consistency condition (|36p yields 



Up to now we have considered the simplest case, where 
f{x, t) is a Gaussian random function with a generic cor- 
relation function. For such a generic random force its 
potential function V{x,t), defined by / = —dV/dx ex- 
ecutes a random walk, exhibiting increasing excursions 
as \x\ — > oo. This generic case corresponds to the case 
where ^ = 1 in (|26p . For general values of C, the long-time 
diffusion coefficient and the law of anomalous short-time 
diffusion in regime II were calculated by [l3| , where it was 
shown that a force derived from a statistically stationary 
potential corresponds to C = 3. 

In this section we briefly discuss the conditions of valid- 
ity for the Fokker-Planck equation for general values 
of (. The discussion follows that at the end of Sec. IIVI 

The three conditions to be satisfied are now: 



• Amplitude condition. From (p7)) we find that 
the steady-state variance of momentum is (p^) ~ 
{poD(:;/^)~'^^^'^''~^^\ The amphtude condition is 
therefore: 



Ap a am^ 
p p p 



x(-l) «1 (49) 



• Frequency condition. In the steady state we require 

UJ \1/(2C+1) 



< 1 



2m^a /cjn2/(2C+i) 



This is the same as (HQ 



(51) 



For C = 1 these three conditions are equivalent and 
correspond to the condition uj'^/x 1 derived in Sec. 
rvl But for other values of C, the conditions are no longer 
equivalent. Condition (|5ip can be written as 



3/2-C 



(52) 



This corresponds to a line with slope — (C— f )~^ in Fig-IH 
For C < 3/2, this condition is satisfied everywhere in re- 
gion II, and therefore does not pose an additional con- 
straint on the validity of the Fokker-Planck equation (|23p . 
For ( > 3/2, by contrast, is not satisfied everywhere 
in regime II; a new non-ergodic region appears in the 
phase diagram above the dividing line ujx~^^^~^'' =const. 
We remark that the condition ([5D|) is always satisfied in 
regime II. 
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